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We study structural and thermophysical properties of a one-dimensional classical fluid made of 
penetrable spheres interacting via an attractive square-well potential. Penetrability of the spheres 
is enforced by reducing from infinite to finite the repulsive energy barrier in the pair potentials 
As a consequence, an exact analytical solution is lacking even in one dimension. Building upon 
previous exact analytical work in the low-density limit [Santos et al, Phys. Rev. E 77, 051206 
(2008)], we propose an approximate theory valid at any density and in the low-penetrable regime. 
By comparison with specialized Monte Carlo simulations and integral equation theories, we assess 
the regime of validity of the theory. We investigate the degree of inconsistency among the various 
routes to thermodynamics and explore the possibility of a fluid-fluid transition. Finally we locate 
the dependence of the Fisher- Widom line on the degree of penetrability. Our results constitute the 
first systematic study of penetrable spheres with attractions as a prototype model for soft systems. 



I. INTRODUCTION 



Hard spheres constitute a paradigmatic system for many simple and complex fluids. Steric stabilized colloids, for 
instance, are suspensions made of colloidal particles coated by short linear polymers suspended in a microscopic solvent 
fluid. For sufficiently high temperature and/or in the presence of a good solvent, those dressed colloids effectively 
interact as hard spheres. 1 
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On the other hand, a number of soft colloidal systems are always penetrable at least to a certain extent 2 . Notable 
examples include for instance star-shaped 3 or branchcd-shapcd 4 polymers where each macromolecue can be roughly 
regarded as a sphere of a given radius (the radius of gyration), but two particles can clearly interpenetrate to a 
substantially smaller distance. 

A necessary (but not sufficient) condition for a one-dimensional fluid to be a nearest-neighbor fluid is to be a 
hard-core fluid, i.e., a fluid made of particles which cannot penetrate one another due to the existence of an infinite 
repulsive potential barrier in the pair potential <f>(r). Nearest-neighbor fluids admit an analytical exact statistical- 
mechanical solution: 5 the partition function, equation of state, and correlation functions of any order can be calculated 
analytically from the knowledge of the pair potential. This is no longer the case for non-neighbor fluids. 6 

Penetrable spheres (PS) 7,8 can be reckoned as the simplest representation of soft colloids where the range of 
penetrability can be tuned from zero (hard spheres) to infinity (ideal gas). Both limits are amenable to an exact 
analytical treatment, but the intermediate case is not. 

When an attractive, short-range, square well (SW) is added to PS, one obtains the so-called penetrable-square- 
well (PSW) fluid. 10 On the one hand, this enriches the model so that it can also account for short-range attractive 
interactions which arc ubiquitous in such systems. On the other hand, it also complicates the treatment due to 
possible Ruche instabilities associated with the lack of a well defined thermodynamic limit. 11,12 As the width of 
the well vanishes with a constant area under the well the PSW model reduces to what we denote 10 as the sticky- 
penetrable-sphere (SPS) model. This model was found to be thermodynamically unstable 10 due to the divergence of 
the fourth virial coefficient. In fact, SPS model violates the (sufficient) condition for stability (see Appendix A in Ref. 
10). 

We emphasize that various classes of penetrable systems have appeared in the literature with rather different 
meanings. The Widom-Rowlinson model of non-additive hard-sphere mixtures, 13 for instance, is not associated with 
a well defined pair potential as in the case of the present study. Likewise, the Rikvold-Stell-Torquato "permeable 
sphere" model 14,15 is defined through a condition on correlation function which is not equivalent to a constant repulsive 
potential inside the core region. On the contrary, our PSW model belongs to the same class of bounded potentials 
as the Gaussian-core models originally proposed by Stillnger et al. 16 in the late 1970s and exploited more recently by 
the Diisscldorf 8,9 and the Cambridge 17 groups. 

In a previous paper, 10 we have introduced the PSW fluid model and discussed the conditions under which the model 
is Ruelle stable. In addition, we have also derived an exact low-density expansion up to second order in the radial 
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pair distribution function (corresponding to the fourth order in the virial coefficient) which was shown to compete 
with standard integral equation approximations such as Percus-Yevick (PY) and hypcrnctted chain (HNC) over a 
wide region of the density-temperature phase diagram. These exact results, however, fail to reproduce the correct 
behavior when the concentration is large, due to their low-density character. 

The aim of the present paper is to extend the analysis to these more demanding conditions, by using an approxima- 
tion already successfully exploited in the PS case. In this case it has been argued that the exact analytical solution 
stemming from corresponding hard-sphere particles can be efficiently exploited to implement a low-penetrability 
approximate solution (LPA, called LTA in Ref. 7). The basic idea behind the method is that for sufficiently low 
penetrability, the functional form of the equations derived in the impenetrable case can be smoothly adapted to the 
penetrable case by "healing" a few crucial aspects of the original solution. Building upon this idea, we here show that 
this methodology can also be applied to the PSW case by starting from the corresponding impenetrable counterpart 
(i.e., the SW potential). 

We discuss the soundness of this approximation in various ways. First by comparing the LPA low-density results 
against the exact low-density expansion which was computed in Ref. 10. Secondly, by comparing with specialized 
Monte Carlo (MC) simulations and standard integral equations (notably PY and HNC). We show how LPA properly 
describes a significant part of the phase diagram with a performance comparable with integral equations at a semi- 
analytical level. 

The introduction of an attractive part in the PS potential opens the route to some interesting questions that we 
also address in the present paper. First of all, we question the existence of a fluid-fluid phase separation in addition to 
the fluid-solid transition, by limiting our analysis within the range of applicability of LPA, that is, we avoid densities 
so high that a substantial interpenetration among particles is expected. 

Within the same LPA, we also investigate modifications on the Fisher-Widom line, marking the transition from os- 
cillatory to exponential decay regimes for correlation functions, that is known to exist even in the SW one-dimensional 
fluid. 18 We find an increase of the exponential decay region and we address the physical motivations behind this. 

The structure of the paper is as follows: we define the PSW model in Sec. II. In Sec. Ill we briefly recall the well 
known general scheme allowing for the exact analytical solution of the class of nearest-neighbor one-dimensional fluids. 
We then construct the LPA in Sec. IV and show how this reduces to its counterpart within the PS limit 7 and assess 
its performance in comparison with known exact results within the low-density limit. 10 Sections V and VI contain a 
discussion on the Fisher-Widom line and on the routes to thermodynamics, as predicted by the LPA, respectively. 



The regions in the density-temperature diagram where the LPA is only slightly thcrmodynamically inconsistent (and 
thus expected to be reliable) are discussed in Section VII, where also an improved version of the approximation is 
proposed. Section VIII includes a very brief description on the numerical methods (MC simulations and integral 
equations) discussed in the present model. These numerical results are presented and compared with LPA theory in 
Section IX. The paper ends with some concluding remarks in Sec. X. 

II. THE PENETRABLE-SQUARE- WELL (PSW) MODEL 

The PSW fluid is defined through the following pair potential 10 (see Fig. 1, top panel) 

:r, r < a, 

4>{r) = { _ £q; a <r <(j + A, (!) 
0, r>a + A, 

where e r and e a are two positive constants accounting for the repulsive and attractive parts of the potential, respec- 
tively. The corresponding Mayer function f(r) = e~' 3l ^ r ) — 1 (where (3 = 1/ksT is the inverse temperature parameter) 
reads 

/ (r) = Tr/ns (r) + 7« [6 (r - a) - Q (r - a - A)] , (2) 

where j r = 1 — e~ /3er is the parameter measuring the degree of penetrability varying between (free penetrability) 
and 1 (impenetrability) and -f a = e t3ta — 1 > plays a similar role for the attractive part. Here JhsM — 0( r — f) — 1 
is the Mayer function for the hard-sphere model which can then be recovered in the limit j r — > 1 and either 7 a — ► 
or A — * 0. 0(r) is the usual step function equal to 1 for r > and otherwise. It also proves convenient to introduce 
the ratio 7 = j a /"fr, which is a measure of the relative depth of the attractive well. 

The above potential reduces to the corresponding PS and SW potentials in the limits e a — ► (or A — > 0) and 
e r — > 00, respectively. Other interesting limiting cases have already been detailed in Ref. 10. 

Consider a SW fluid in one dimension: different particles can be assigned an increasing coordinate on the axis line 
and the only possible configurations arc those indicated with A or B in Fig. 1 (middle panel) , where cither the centers 
of two different spheres are separated a distance greater than the attractive SW range and behave as hard spheres 
(A) or they are sufficiently close to attract each other (B). PSW spheres, on the other hand, can interpenetrate with 
some energy cost so they also display configurations such as, for instance, (C) or (D) in Fig. 1 (bottom panel). PSW 
fluids are then effectively a many-body problem and, as such, not amenable to an analytical solution. 
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In the present paper, our analysis will be limited to the case e r > 2e a where a well defined thermodynamic limit is 
ensured. 10 



III. GENERAL RECIPE FOR NEAREST-NEIGHBOR INTERACTIONS 

In this section we provide a synopsis of the main steps required by the analytical solution of any nearest-neighbor 
fluid. 5,19,20 This will be used in next section to introduce a motivated approximate solution in a particular limit. 

• From the Boltzman factor e _ ^ (r <' compute its Laplace transform 

/•OO 

fi(s) = / dre- sr e-^ {r) . (3) 
Jo 

• The equation of state is given by 

/%>=-, (4) 
a 



(5) 



where p is the pressure and the parameter £ is the solution of the equation 

where p is the density and fi'(s) = d£l(s)/ds. This provides all thermodynamics. 

• The radial distribution function (RDF) can be obtained from 

GOO = ~~ " (B t^> , (6) 

p n (f/tr) - n (s + e/tr) 

which is the Laplace transform of the RDF <?(r). 

This is sufficient to compute both thermodynamics and structural properties of any one-dimensional system with 
nearest- neighbor interactions. 

At odds to this class of problems, penetrable spheres do not possess any analytical solution even in one dimension. 
This is because it is not possible to convolute appropriate Laplace transform along a one-dimensional axis, which is 
the essential feature rendering the short-range one-dimensional models solvable. In turn this is due to the existence 
of multiple "blobs" formed by interpenetrating spheres so that it is no longer possible to "order" them along a line 
in such a way that they do not cross each other, a key point to the existence of the analytical solution (see Fig. 1, 
middle panel). 

Because of this, we now turn our attention to a motivated approximation which amounts to assume a slight decrease 
from an infinite repulsive barrier, an approximation which will be denoted as low-penetrability. 
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IV. THE LOW-PENETRABILITY APPROXIMATION (LPA) 
A. Construction of the approximation 

In Ref. 10 we have followed the philosophy of considering a low-density expansion to provide exact analytical results 
valid up to second order in the RDF g(r) and up to fourth order in the virial expansion. This is a very useful exact 
limit case to test approximate theories and numerical simulations, but it has the considerable disadvantage of being 
limited to very low densities. We now consider a different approach where density can in principle be arbitrarily large 
but we assume low penetrability among different spheres, patterned after a similar idea already used in the PS case. 7 

For notational simplicity, in the following, lengths will be measured in units of a (so that a = 1) and we introduce 
A = 1 + A/cr as a dimcnsionlcss measure of the external well boundary. 

The Laplace transform of the Boltzmann factor e^ 1 "-' for the PSW model is 

fi(«) = i-^ + ^[(l + 7)e -- 7e ^]. (7) 

The PSW fluid is not a nearest-neighbor fluid, as remarked, but it reduces to the nearest-neighbor SW fluid as 
7r — > 1 (and A < 1). In this limit, it is natural to use the recipe given in Sec. Ill for the SW fluid, to derive an 
approximate equation of state and an approximate g(r) from Eqs. (4), (5), and (6), respectively. This, however, 
must be exercised with care as important general properties of any model, such as for instance the continuity of the 
cavity function y(r) = g(r)e^^ r \ are typically lost by this brute force procedure. The driving idea behind this simple 
low-penetrability approximation (LPA) is then to keep the general features of the original SW solution and enforce 
some specific modifications guided by the accounting of increasingly important constraints. 21 
Our LPA implementation amounts to replacing Eq. (6) with 

G(s) = ^ + , (8) 

pO (C)-O (s + C) 

where flo(s) is (formally) the Laplace transform of the Bolzmann factor of the SW model which can be obtained from 
f2(s) by discarding the first term on the right-hand-side of Eq. (7), i.e., 

h (s) = 2l[(i + 7 )e-- 7e -*»]. (9) 
s 

This simple choice can be shown to be fully equivalent to keeping Eq. (6) but with a replacement e~ s — * e~( s ~ a ^ in 
Eq. (7), where the free parameter a is fixed by the continuity condition of the cavity function y(r) at the hard-core 
discontinuity r = 1. This is known to be the most important feature to obtain a correct representation in integral 
equation theories of SW fluids, both from the analytical and the numerical viewpoint. 22-24 



We note that, unlike the SW counterpart, C 7^ Pp. It is a transcendental function of (3 and p which can be obtained 
by ensuring the correct behavior of g(r) — > 1 as r — > 00 or, equivalently, sG(s) — > 1 as s — > 0. From Eq. (8), this gives 

" (0 -C k + ]-\^ do) 



where in the second equality we have introduced the following quantities: 

« - 7^, (11) 
h = -^TtV" (12) 

7r(l + 7) 

For given values of the potential parameters (A, e r , and e a ) and for given values of the inverse temperature /3 and 
the auxiliary parameter £, the quantities g and /i are obtained from Eqs. (11) and (12) and inserted into Eq. (10) to 
determine the density p. The impenetrable SW potential corresponds to the limit h — > 0. 
In order to compute the RDF g(r) we first compute explicitly G(s) from Eqs. (6) and (7) 



h + e- s (1 - qe~ sA ) 
p (1 - q){l + s/()~e- s {l-qe- sA ) ' 



G(s) = ^ _ w /, ^ ^ (13) 



Upon expanding the denominator in Eq. (13) in powers of (1 — qe )(1 + s/C), and inverting the Laplace transform 
term by term one gets 



, j. 00 n / 

n=l fe=0 

where 



W M = T^e-f + E E J (-^)^„ (r - n - kA) 6 (r - n - kA) , (14) 



c 



1-9 



n r ^"-1 



e-< r . (15) 



(n-1)! 

We anticipate that the LPA does not capture correctly the r < A trend at high densities, while it works well for 
r > A. The reason for this can be traced back to the failure of the LPA to account for the discontinuous slope of the 
cavity function y(r) at r = A. Moreover, the appproximate y(r) turns out to be discontinuous rather than continuous 
at r = A, as detailed in Appendix A. These deficiencies can be accounted for step by step at the price of an increase 
in the complexity of the approximation, and are a consequence of the phcnomcnological nature of the LPA. This will 
be further discussed in Sec. VII. 

As already remarked, the PSW model reduces in the appropriate limit to the penetrable analogue of Baxter's sticky 
hard spheres, denoted as SPS in Ref. 10. This is further elaborated in Appendix B, where it is also discussed the 
LPA of the SPS model. We have explicitly checked this is indeed the limit for PSW in the limit of very narrow and 
very deep well. On the other hand, we have also found (see Appendix B) that this model is also thermodynamically 
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unstable as it violates the stability criterion e r > 2e a , and hence it will not be further discussed in the remaining of 
this paper. 

B. The penetrable- rod limit 

Here we show that either in the limit e a — > (which implies 7^0) or. alternatively, in the limit A — > 0, the LPA 
that we just found for the PSW model reduces to the corresponding one proposed in Ref. 7 for the PS model. 

Taking the limit 7 — > in Eq. (7) one finds Eq. (2.53) of Ref. 7. Moreover q — > and h — > (77 1 — l)eS and so Eq. 
(10) reduces to p = [1+ (77 1 — l)e^]/(l + C _1 ) which can be rewritten as (£ — C) e_ V ilr 1 ~ 1) = C with £ = (0(1 + £), 
which coincides with Eq. (4.4) of Ref. 7 where our £ replaces their 

It is straightforward to check that the same expressions for f2(s) and for p in terms of £ and 7 r are obtained in the 
alternative limit A — > 0. Hence LPA for PS is fully recovered. 



C. Comparison with exact low-density expansion 

It proves interesting to compare the LPA to order p with the exact results derived in Ref. 10 based on a low-density 
expansion, in order to assess the ability of LPA to reproduce low-density results. 

The general expansion of g(r) in powers of the density p has the following structure 25 



9 (r) = 90 (r) + gi (r) p + 



(16) 



The exact results for go(r) and g\{r) have been derived in Ref. 10: 

l-7r, r < 1, 

^cxact f„\ ) 

9o (r) - <i + 77r; 1 < r < 1 + A, 



(17) 



1, 



r > 1 + A, 



(1 - 7 r) [2 (1 + 7 2 A) - r (1 + 2 7 + 2 7 2 )] , < r < A, 
(1 - 7 r ) (2 - 2 7 A - r) , A < r < 1, 

(l+77r)(2-27A-r), 1< r < 1 + A, 



ff r act (r) - 7 2 



2 - 27A - r, 

7(2 + 7) (r-2)-2 7 A, 

(2 + 2A-r) 7 2 , 

0, 



1 + A < r < 2, 



2 < r < 2 + A, 



2 + A<r<2 + 2A, 



2 + 2A < r. 



(18) 



In order to compare <?o xact ( r ) an d 3i xact ( r ) with LPA results, we expand £ as derived from Eq. (10) to lowest order 
in density, C = CoP + CiP 2 + ^>(p 3 ), and plug the results into Eqs. (14) and (15). This yields Eq. (16), where the 
coefficients <?o( r ) an d gi(r) are computed within the LPA. Whereas go(r) = <?o xact ( r )' 9i( r ) ^ s found to differ from the 
exact result <?J xact (r). Analytical expressions for (0, £1, and gi(r) can be found in Appendix C. 

Having done this, one can estimate the difference in the cavity function between LPA and exact results to order p, 
which reads (see Appendix C) 



yr C \r)- yi {r) = 



Ci+L>iA + £i(A-l) + Fi(r-A), < r < A, 
Ci+L>iA + £i(r-l), A<r<l, 



(19) 



where 



C 1 + D 1 (l + A-r), 



d = 77 2(i_ 7r )J±T_A, 

1 + 77r 

Di = 77r— , 

1 + 77r 
El = 7r(l ~7r), 

Fi = 7r[l-7r-277 r (l + 7 )]. 



1 < r < 1 + A, 



(20) 

(21) 
(22) 
(23) 



The right-hand side of Eq. (19) preserves the continuity of yi(r) at r = A and r = 1, but imposes the continuity of 
yi(r) at r = 1 + A and that of y[(r) at r = 1 and r = 1 + A, as well as the discontinuity of the exact y[ (r) at r = A. 
The latter discontinuity is, according to Eqs. (17) and (18), 

y'(A+)-y'(A-) 



lim 



= 277^(1+7). 



(24) 
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v. 



THE FISHER- WIDOM LINE 



In a remarkable piece of work, 18 Fisher and Widom argued that the asymptotic decay of the correlation functions 
is determined by the nature of the poles = Sj(/3, p) (i — 1,2,3,.. .), with largest real part, of the Laplace transform 
G(s) of the RDF. This asymptotic decay can be of two different types: oscillatory at high densities and/or high 
pressures, and monotonic for low densities and/or pressures. The latter regime can exist only in the presence of 
competing effects in the potential function, so it cannot exist for purely repulsive short-range potentials, such as HS 
and PS potentials. 

In particular, rather general arguments 26 suggest a behavior 



where we have specialized to one-dimensional systems and the sum runs over the discrete sets of poles Sj, A4 being (in 
general complex) amplitudes. The asymptotic behavior of g{r) is dominated by the pole s\ having the least negative 
real part (to ensure stability of the liquid). If s\ is complex, its conjugate S2 — s\ must also be included in the 
asymptotic behavior. 

Fisher and Widom derived the line — henceforth denoted as Fisher- Widom (FW) line — where this transition takes 
place, both in the pressure versus temperature and in the density versus temperature diagrams, for the one-dimensional 
SW potential. On crossing this line, one finds a sharp transition in the character of the RDF g(r) — 1: For any fixed 
temperature in the p-T plane, g(r) has an oscillatory character above the FW line and an exponential decay below 
it. The transition is a signature of local ordering without any singularities in thermodynamical quantities as there is 
no phase transition in the one-dimensional SW fluid. In three-dimensions, the FW line precedes the coexistence line 
when lowering the pressure at a fixed temperature. This has been numerically observed for various fluids including 
SW, 26 Lennard- Jones, 27-29 and other softer potentials. 30 

In view of the possibility for PSW to display fluid-fluid and fluid-solid phase transitions in spite of their one- 
dimensional character, it is interesting to wonder what happens to the FW line in the transition from SW to PSW. 
We now analyze this in the framework of the LPA. 

The poles of G(s) (different from s = 0) can be read off from Eq. (8): 



As we are here interested in the pole with the negative real part closest to the origin we set s = —x 7^ as the real 




(25) 



n (s + () = n (0 . 



(26) 
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root of Eq. (26), 

n Q (C - x) = n (c) , (27) 

and s = —x' ± iy as the complex root with the least negative real part, i.e., 

Ren (C-x'±iy) = fio(C) , (28) 
lmh (( - x' ± iy) = . (29) 

The pole s\ determining the asymptotic behavior is cither si = — x (monotonic decay) if x < x' or si = —x' ± y 
(oscillatory decay) if x > x'. The FW transition takes place when x = x' . 
Equation (27) yields the condition 

e- s (l-qe- sA ) = (1 - q) (l + , (30) 

where q is given by Eq. (11). Quite interestingly, as the parameter h is missing, this equation formally coincides with 
its SW counterpart, originally studied by Fisher and Widom [see Eq. (3.6) in Ref. 18]. We can rewrite Eqs. (27)-(29) 
as follows: 



e x(l- q e^) = (l-q)(l--), (31) 
e x ' (cosy - qe x ' A cos\y) = (1 - g) - , (32) 
e x (sin y - qe xA sm\y) = -(1 - ? )|. (33) 

At the FW transition (x = x'), Eqs. (31)-(33) form a set of three coupled equations whose solution yields x, y, and 
£ as functions of q. Use of Eq. (10) then gives the line in the p-T plane. 
It proves convenient to eliminate £ from Eqs. (31) and (32) to obtain 

x = ~r^—, Ti (34) 

so that from Eqs. (31) and (33) we can now get 

cosy-cosAy 

Q = x — y — . (35) 

sin y — sm Xy + sin y A 



When Eq. (34) and (35) are inserted into Eq. (31) we get 



y — y 

sin y — sin Xy + sin y A (cos y — cos Xy) — — e x (l — q)— (1 — cos Xy) , (36) 

x x 
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where x(q, y) is given by Eq. (34) so that this is a transcendental equation in y(q). Once y{q) is known from Eq. (36), 
Eqs. (34) and (35) provide x(q) and £(<?), respectively. The parameter 7(5) is obtained by inverting Eq. (11), 

= e -cjA_ g - ( 3? ) 

and the inverse temperature (3{q) is obtained from 

= 1 - e -«.K (38) 

on using the definitions of 7, 7 r , and 7 a . 

Finally, Eqs. (10) and (12) provide and the combination of j3(q) and p(q) yields the FW line in the p-T plane. 
In order to have it in the p-T plane one needs to get before the equation of state and the result will depend on the 
chosen route (virial, compressibility, or energy). This is discussed in the following section. 



VI. EQUATION OF STATE 



As PSW is not an exactly solvable model, thermodynamics will in general depend upon the followed route, so we are 
going to check the three standard routes (virial, compressibility, and energy) for the compressibility factor Z = (3p/ p, 
as predicted by the LPA. 

The virial route is defined by 



Z = / dr ry (r) e~ mr V (r) 

Jo 



which, using standard manipulations, yields 



Z = l + p 7r [(l + 7)y(l)-7Ay(A)]. 



(39) 



(40) 



As y(A ) ^ 2/(A + ) within LPA (sec Appendix A), y(X) = | [y(X ) + y(X + )] is assumed. Thus, using Eqs. (A4) and 
(A8) we get 



Z = 1 + 



1-5 



1 - Xq 



l + 7r(l + 7) 



1 + 7 7r /2 h(A 



(41) 



I + 77,. l-<7. 

It is easy to check using Eqs. (10) and (41) that in the case of the SW model (h = 0) one recovers the expected result 

z = a p. 

Next we consider the compressibility route: 



= l + 2plim [g(s)-s~ 1 



(42) 



Using Eqs. (10) and (13) the last term of Eq. (42) can be written as 

'd P \ h'(0 



2p Iim 

s->0 



G (s) - s" 



Introducing the quantity 



Eq. (42) becomes 



xin = i n(,(C)-n'(C) 
^ ~ (dp/dOp Q' (0 : 



and using the definition of x we nn d 



ft 



(dp) _i 



Therefore the compressibility route yields 



f c dC 

MO = 



In the SW limit one clearly has X(£) — and (3p — £, as it should. 

The energy route is by far the most delicate. We start from the internal energy per particle 

1 f°° 
u = — + p I dr 4>(r)g(r) 



2/3 

= ~^j + er pf o dr 9(r)-e a pJ^ drg{r). 
Equation (Al) provides the necessary result for g(r) in the interval < r < A, so that 



1 h , _> N 



l-q V '\ 1-q J (1-9) 2 

In order to obtain ftp from u we exploit the following thermodynamic relation 



and the identity 



to obtain 



dp) (3 \dUf3\dP/p 



fd/3p\ p 2 /du\ 
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Once again one can check that Eq. (52) is satisfied by the SW result (3p = (. 

The right-hand side of Eq. (52) is a function of (3 and p, which we denote as R((3, p), as £ is itself a function of the 
same variables through Eq. (10). Thus, Eq. (52) gives 



We are now in the position to draw a qualitative phase diagram in the pa-ksT / e a plane indicating the boundary 
where the LPA can be approximately regarded to be reliable. Of course, a definite reliability test is only possible after 
comparison with computer simulation results but before that we can use the internal consistency among the three 
thermodynamic routes as a reliability criterion. 

In general, it turns out that thermodynamic inconsistency increases as the temperature and the density increase. 
To characterize this, let us define a density pu m (T) such that the largest relative deviation among the three routes is 
smaller than 5% if p < piim(T). Therefore, all the points in the temperature-density plane with p < pn m (T) represent 
states where the LPA is only weakly inconsistent. This boundary line is shown in Fig. 2 for three representative cases 
of the pair (e r /e a and A). We observe that the region where the LPA is thcrmodynamically consistent shrinks as e r /e a 
decreases and/or A increases. In any case, it is noteworthy that if the density is smaller than a certain value (which 
of course depends on e r /e a and A), the LPA remains thcrmodynamically consistent even for high temperatures. 

The above reliability criterion is based on thermodynamics and thus it is a global one. On the other hand, we 
know that the LPA has some local shortcomings, such as an artificial discontinuity of the cavity function at the point 
r = 1 + A, as shown in Appendix A. Moreover, it docs not predict the discontinuity of the slope of the RDF at r = A, 
already present by the exact result to first order in density, as indicated by Eq. (24). 

As anticipated in Sec. IV, we can extend the validity of the LPA by a suitable modification of the cavity function 
y(r) in order to ensure a correct behavior both within the core region and at the well-edge discontinuity. We outline 
a possible approach to this issue in the remainder of this Section. 

Inspired by the comparison with exact low-density results as given in Section IV C, we modify the LPA (mLPA) 




(53) 



where /3, 



max 



is a conveniently chosen high value. 
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by adding linear terms in the region < r < 1 + A, following a form based on that of Eq. (19), namely 



mLPA 



(r) = g(r) + -{ 



(1 - 7 r )[C* + DA + E(A — 1) +F(r — A)], < r < A, 

(l-7r)[C + DA + E(r-l)], A<r<l, ( 54 ) 

(l+7 7r )[C + L>(l + A-r)], l<r< 1 + A, 

where g(r) is the LPA radial distribution function as given Eq. (14). The parameters C, D, and E can be determined 
by imposing the continuity of y(r) at r = 1 + A and of y'(r) at r = 1 + A and r = 1, respectively. They are given by 

7r(l + 7) ^C 2 



c = 

D = 

E = 



i + 77r (i-g) 2 
i - 7r fe g c 2 
i + 77r(i-g) 2 
l /iC 2 



A. 

-cc, 

L>. 



(55) 
(56) 
(57) 



l+ 77r (1 -g) 2 

The addition of the coefficient F is motivated by the exact results to first order in density, Eq. (18), showing that, 
as recalled above, g(r) exhibits a change of slope at r — A, a feature not accounted for by the LPA. In order to 
determine the coefficient F we extend the exact low-density condition (24) to finite density. This implies 

1+7 hq( 



F = E-2fr)Y r {l+l) 



C + DA + E(A - 1] 



(58) 



.7(1 -7r) 1 -q 

It is straightforward to check that C = dp 2 + 0(p 3 ), D = D lP 2 + 0(p 3 ), E = E lP 2 + 0(p 3 ), F = F lP 2 + 0(p 3 ), 
where Ci, D\, E\, and Fi are given by Eqs. (20)-(23). Therefore, the mLPA is exact to first order in density. 

The discussed modification of LPA then takes care of the continuity of the cavity function y(r) at both interaction 
discontinuities r = 1 (already accounted for within LPA) and r = A (where the original LPA fails to provide continuity), 
and it correctly matches the exact results for g(r) up to first order in density. A similar modification of the SPS 
model, discussed in Appendix A, would heal the discontinuity appearing in the corresponding LPA values y (1 + ) 7^ 
y SPS (l~), which is a consequence of the combined effects of the LPA discontinuity y(X + ) 7^ y(^ ) and the sticky 
limit. This would provide an expression (not reported here) which is this sticky limit of Eq. (54). 



VIII. MONTE CARLO SIMULATIONS AND INTEGRAL EQUATION THEORY 



In order to assess the reliability of the LPA, we will compare in Sec. IX with specialized MC simulations. In 
addition, prompted by the results of Ref. 10, we will also compare LPA with standard integral equation theories, such 
as PY and HNC. 25 
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A. 



Monte Carlo simulations 



We have employed the conventional Monte Carlo simulation on an NVT ensemble with periodic boundary conditions 
which in one dimension means that the system is treated as a ring. N = 5 x 10 penetrable-rod particles were displaced 
according to the Metropolis algorithm to create an initial sample of configurations. Following the equilibration stage, 
each run is divided into 20 basic simulation blocks, in which 10 5 measurements are performed to collect correlation 
functions data. 100 trial moves per particle are implemented between each measurement, so that 10 13 equilibrium 
configurations arc generated in each run. 

In order to speed up the simulation process the particles are labeled such that they create a consecutive sequence 
in a clockwise order. Calculation of a potential of a particle i in a given configuration then reduces to a searching for 
the highest label j > i and the lowest label k < i associated with the particles still interacting with the particle i. 
In contrast to the case of impenetrable spheres in one dimension the order of particles changes so that a relabeling 
must be undertaken after each shift of a particle. Obviously, at higher temperatures the number of penetration can 
be high, which makes the calculations more demanding compared to hard body systems. 

There are in general two routes for the evaluation of the pressure. Determination of the pressure using a mechanical 
(virial) route relies on an ensemble average of a virial, i.e. a quantity involving the forces acting on all the particles. 
Alternatively, a thermodynamic expression relates pressure to the volume derivative of the free energy and is imple- 
mented by calculating the free energy change associated with small virtual change of volume. However, for systems 
with discontinuous interaction both mechanical and thermodynamic approaches become identical. Specifically, for the 
PSW fluid model both approaches reduce on a calculation of distribution function at r = 1 and r = A [sec Eq. (40)]. 



The presence of penetrability does not pose any particular difficulties to standard integral equation theories. As a 
matter of fact these have been already employed in the PS case 7 and in the PSW case 10 within standard approximations 
where the one-dimensional Ornstein-Zcrnike equation 



B. Integral equations 




(59) 



is associated with a PY closure 



c(r) = f{r)y{r) 



(60) 
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or with an HNC closure 

c (r) = / (r) y(r)+y(r) — l — In y (r) . (61) 

We have solved the PY and HNC integral equations using a Zerah's algorithm 32 with up to 2 12 grid points depending 
on the considered state point. 

IX. RESULTS WITHIN THE LPA 

In this Section we compare numerical results stemming from the LPA with MC simulations and integral equation 
theories (PY and HNC) for both RDF (where we will consider the improved mLPA) and equation of state (at the 
level of the simple LPA). 

A. Results for g(r) 

As a first approach to assess the performance of the LPA, we consider the RDF g(r) for two representative state 
points. 

The well is kept fixed at A/a = 0.5 and temperature is also fixed by the attractive energy scale so that fcsT/e a = 1. 
Figure 3 depicts the behavior of g(r) for a density pa = 0.2 and an energy ratio e r /e a = 5, which is well above 
the stability threshold value e r /e a = 2. 10 The stability threshold is then probed in Fig. 4, whereas a higher density 
pa = 0.8 is tested in Fig. 5 with all other parameters identical to those of Fig. 3. 

In all cases, mLPA results (that only differ from the LPA ones within the interaction range, < r < A) are 
compared with MC simulations and integral equations and follow the expected trend. For low densities (pa = 0.2) 
and low penetrability (e r /e a = 5) mLPA, PY, and HNC all provide very accurate descriptions of MC data with a 
very tiny difference in the well region 1 < r/a < 1.5, where the integral equations predict a slight curvature of g(r), 
while the mLPA confirms the practically linear shape of the simulation data. Moreover, a blow-up of g(r) in the deep 
core region (0 < r < A) shows that the mLPA is very accurate, while the PY and HNC theories underestimate and 
overestimate, respectively, the MC data. The same good performance of the mLPA is also observed for a much larger 
penetrability (e r /e a = 2), provided the density is relatively low (pa = 0.2), as shown in Fig. 4. This is consistent 
with Fig. 2, according to which the density pa = 0.2 lies in the region where the LPA is expected to be accurate for 
any temperature when e r /e a = 2 and A/a = 0.5. As for the integral equations, they are also rather accurate for the 
case considered in Fig. 3, although they still show a slight curvature inside the well and slightly deviate from the MC 
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results for r < A. Differences begin to be relevant at high-density (pa = 0.8), mostly inside the core < r/a < 1 
and in the contact values r = a + . Again, this agrees with Fig. 2, which shows that the state (per, fcsT/e a ) = (0.8, 1) 
is practically on the boundary line corresponding to e r /e a = 5 and A/a = 0.5. In any case, Fig. 5 shows that the 
best general agreement with the MC results is presented by the mLPA, followed by the HNC theory, which, however, 
predicts reasonably well the peaks of g(r), but not the minima. 

We have explicitly checked (not shown) that for smaller values of the well width A, PSW results increasingly tend 
to the SPS counterpart, as anticipated. 

B. Results for equation of state 

Next we turn to the analysis of thermodynamics within LPA. As anticipated (see Section VI), the lack of an exact 
solution gives rise to thermodynamical inconsistencies where compressibility, virial, and energy routes all give rise to 
different results. The consistency degree among different routes is a (partial) signature of the LPA performance, as 
discussed in Sec. VII. In Fig. 6 we report the behavior of 0p as a function of the reduced density pa. Once again, wc 
fix the width of the well A = 0.5<r and the energy ratio e r /e a = 5, and consider two different temperatures fcgT ' /e a = 1 
(top panel) and ksT/e a = 5 (bottom panel). In the former case different routes give practically indistinguishable 
results up to pa ~ 0.8, whereas in the latter a difference is clearly visible at densities higher than pa ~ 0.6 with energy, 
virial, and compressibility routes having decreasing (3p for identical values of pa. Similar results are observed at the 
stability edge e r /e a = 2, as shown in Fig. 7. We remark that higher temperatures effectively correspond to higher 
penetrability, as particles have relatively more attractive energies, as compared to the positive repulsive barrier, and 
hence they can compenctrate more. Therefore pressure differences among different thermodynamical routes can be 
reckoned as a rough measure of the breakdown of LPA. On the other hand, consistency among different routes does 
not necessarily means "exact" results, as they can all converge to the incorrect value. 

A comparison with MC numerical simulations is therefore also included in Figs. 6 and 7. Somewhat surprisingly, 
this suggests that the virial route is the closest to the true value for the pressure, with both compressibility and energy 
routes always lying on the opposite side with the latter being the farthest from the MC results. 

In order to compare with LPA, we have carefully scanned a wide range of temperatures and densities within the 
region < pa < 1 where LPA provides consistent thermodynamics as remarked. Within this region we found no 
signature of fluid-fluid transition line as expected. Our preliminary numerical results for higher densities, where strong 
overlapping among different particles is enforced, provide a clear evidence of phase separation. As the main emphasis 
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of the present paper is on analytical approximations, this point will be discussed in some detail elsewhere. 

C. Results for Fisher- Widom line 

Let us follow the recipe given in Sec. V to locate the FW line. In Fig. 8 we report the quantities pa/e a and pa as a 
function of ksT/ta for A = a and decreasing values of the ratio e r /e a . The case e r /e a —* oo is the one addressed in 
the original FW work on the one-dimensional SW fluid. 18 We remind that above the FW line, g(r) — 1 has oscillatory 
behavior, whereas it is exponentially decaying below it, and it is located in the homogeneous fluid region of the phase 
diagram, above the critical temperature if phase separation is present. 

As the repulsive barrier becomes finite, the region of monotonic behavior increases for large ksT/ea whereas it 
remains essentially unchanged for lower temperatures. This is not surprising as penetrability (i.e., finite repulsive 
barrier) favors the onset of a critical region. Somewhat more surprising is the fact that this happens in the high- 
rather than in the low-temperature region. A similar feature is also appearing in the p-T plane (see bottom panel). 
In order to test the effect of different width values, we have repeated the same calculation for A = 0.5cr. Results 
are presented in Fig. 9 and are in agreement (in the limit e r /e a oo) with results for the one-dimensional SW fluid 
presented in Ref. 33 for a hard-core to well- width ratio equal to 2 (see Fig. 1 in Ref. 33). For this well width the 
influence of the ratio e r /e a on the FW line is much less important. 

Although we have been unable to find a simple physical explanation for this behavior, we remark that the sensitivity 
of the FW line to the barrier height occurs as the density decreases. Consider for instance the density pa = 0.1 for 
models with A = a. In the SW case (e r /e a — * oo) the decay of the RDF changes from monotonic to oscillatory as 
one increases the temperature and crosses the value kBT/e a ~ 2.2. In the case of the PSW model with e r /e a = 5, 
according to the LPA, the transition takes place at fcsT/e a ~ 2.8. If the density is sufficiently low (pa < 0.076 for 
fr/fa = 5), the asymptotic decay of g(r) — 1 is monotonic for any temperature, while this effect is absent in the 
impenetrable SW limit. One might argue that this influence of the energy ratio e r /e a on the high-temperature branch 
of the FW line is an artifact of the LPA since the latter approximation is a priori restricted to low temperatures. On 
the other hand, this high-temperature branch also corresponds to low densities, counterbalancing the penetrability 
effect and making the LPA presumably accurate. As a matter of fact, the FW lines plotted in the top panels of Figs. 
8 and 9 are obtained from the three thermodynamic routes but the three curves are, in each case, indistinguishable 
each oth! er. In other words, the FW lines are well inside the regions in Fig. 2 where the LPA is thermodynamically 
consistent from a practical point of view. 
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X. CONCLUSIONS AND OUTLOOK 

Onc-dimcnsional fluids with nearest-neighbor interactions admit an exact analytical solution for both structural 
and thcrmophysical properties with a well defined protocol. 5,20 Nearest-neighbor interactions, in turn, require a well 
defined hard-core term in the pair-wise potential preventing superpositions and particle exchanges which is the crucial 
ingredient necessary for the exact solution. The absence of the above constraint, on the other hand, allows the 
presence of critical phase transitions, in spite of the one-dimensional character of the system, which are fully absent 
in the hard-core counterparts. 

Effective pair interactions with a soft-repulsive component are well-known features of polymer solutions and colloidal 
suspensions. 1 ' 2 Among many different model potentials, 2 with various degrees of core softness, penetrable spheres (PS) 
stands out for its simplicity. In this model, the infinite repulsive energy is reduced to a finite one, thus introducing 
an effective "temperature" into an otherwise athermal hard-sphere system. This potential model lacks of attractive 
interactions but these can be accounted for in the penetrable square- well (PSW) companion model where an attractive 
short-range square- well is added to the PS model. 10 

At sufficiently low temperatures, thermal energy cannot overcome the repulsive barrier and penetrability is low, 
whereas at high temperatures different particles can interpenetrate to a significant extent. Hence, within this frame- 
work, low- and high-temperature and low- and high-penetrability terminology can be used synonymously. 

In this work we have studied structural and thermodynamic properties of the PSW model. Using a low-penetrability 
approximation (LPA) akin to that discussed for PS, we have considered rather interesting issues specific of the 
presence of attractive interactions (and thus absent in the PS model) such as fluid-fluid phase separation or the 
existence of a Fisher- Widom line. 18 This is a pseudo-transition associated with a clear-cut change, from oscillatory to 
monotonic, in the asymptotic decay properties of the radial distribution function, as transition line is approached, even 
in those cases where the existence of a critical region is prevented by rigorous theorems (e.g. the SW one-dimensional 
fluid). It requires the simultaneous presence of attractive and repulsive energies and hence it cannot exist for the 
simpler PS model 

Our LPA has been devised to reduce to that of PS in the limit of no well. We have assessed its performance by 
comparing it with exact results 10 in the low-density limit and by comparing with MC simulations and PY and HNC 
integral equation theories for larger densities where exact analytical results do not exist. We found that it reproduces 
a significant portion of the T-p parameter space at the level of pair correlation function, the main difference being in 
the penetrability region < r < a. At odds with its square-well counterpart, PSW thermodynamics depends upon 
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the chosen route in view of the inconsistencies introduced by the LPA. We have quantified the inconsistencies among 
virial, compressibility, and energy routes and discussed how they reflect into the computation of the Fisher Widom 
line. In all considered cases, we found a magnification at large temperatures of the monotonic regime region as 
penetrability increases and a much smaller, if any, modification, at lower temperatures. In all cases the FW line is 
found within the region where LPA is expected to be accurate as thermodynamic inconsistencies are small. Within 
the density region < pa < 1, we have found no sign of a fluid-fluid phase separation, although both fluid- fluid and 
fluid-solid transitions are expected at higher densities. 

In the limit of infinitely narrow and deep well, PSW has been shown to reduce to a penetrable version of Baxter 
adhesive model, 3 which violates the stability condition set for a well defined thermodynamic limit. 10 

As the main weaknesses of LPA for the PSW stems mainly from a non-adequate representation of the penetrable 
region < rja < 1, we have then discussed how a simple modification of the radial distribution function in this region 
gives a significant improvement when tested against MC results under rather demanding conditions. 

This paper is part of an on-going effort on PSW outlined in our previous work. Future work will address a 
complementary approximation (the high-penetrability limit) and its matching with the LPA discussed in the present 
paper, so that the entire parameter T-p-p space can be discussed with some comfortable degree of confidence. This 
will resolve some of the subtle points with no conclusive answer left by the present paper. In addition, a detailed 
investigation of the high density region pa > 1 is underway and will be reported elsewhere.. 
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From Eq. (14) wc have that if r < 2. 



Pg(r) = Y~q & 



0, 

^i(r-l), 



< r < 1, 

1 < r < 1 + A, 



ipx{r - 1) - g<V>i(r - 1 - A), 1 + A < r < 2. 
The explicit expressions of ipi(r) is, from Eq. (15), 
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The continuity condition of y(r) at r = 1 is then given by condition 
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which is identically satisfied, so that 
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However, y(r) is discontinuous at r = A = 1 + A: 
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The jump is then given by 



T (\+\ (\-W 7r 1 + 7 hg( 2 A 
P J/(A + )-y(A ) = — -, 

1 + 77 r (1 - q) 2 

and the value used as an estimate of the point is then given by the average of the left and right limits 
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APPENDIX B: THE STICKY-PENETRABLE-SPHERE (SPS) MODEL 



In this Appendix, wc provide a connection with the SPS introduced in Rcf. 10. This is the penetrable analogue 
of Baxter's sticky-hard-spherc (SHS) well known model. 34 The SPS limit can be obtained by considering the limit 
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A — > and e a — > oo so that a = 7A remains finite, hence playing the role of an adhesivity parameter. We then define 
SPS by the Mayer function 10 

/sps (r) = 7r/sHS (r) , (Bl) 

where 

/sHs(r) = / HS W + a5+ (r - a) (B2) 
is the Mayer functions of the SHS potential and 

S + (r) S lim ^)-e{r-a) 

a->0+ a 

The fluid parameters are then the adhesivity coefficient a > 0, the penetrability coefficient j r , and the density p. 

As anticipated, the SPS fluid is thermo dynamically unstable in the sense discussed in Section II. This can be 
seen both because the required limit does not satisfy the sufficient condition for stability e r > 2e a , 10 and directly 
using arguments akin to those used by Stcll 35 to prove the instability of the original Baxter's model 34 in dimensions 
greater than one. Nonetheless it provides an overall consistency testbench to the performance of LPA within the well 
established framework of SHS. 

In the combined limit 7 — > 00 and A — > with a = 7A, Eqs. (7) and (9) become 



nf s {s)= lr (a+-\e- s . (B5) 



Using the first equality in Eq. (10) it follows that 
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/ = l —^e< . (B7) 

It 



We then use the LPA as given in Eq. (8) to find 



;sps,,» _ f/(s + Q + [a + l/(s + Q]e- 



= u;i/fl-[ a+ iW. + oi.-- - (B8) 

whose inverse Laplace transform yields the RDF, 

00 

P5 SPS (r) = £ ^ PS (r - n) (r - n) , (B9) 



where 
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V<o SPS M = —rr^e-^ , (BIO) 



a 



e~ Cr . (Bll) 



« + VC, 

In the impenetrable limit j r — > 1 and / — > 0, Eqs. (B6)-(B11) reduce to the exact one-dimensional SHS 
counterparts, 22,36 ' 37 as they should. 

A word of caution is in order here. Using Eqs. (B9)-(B11), the cavity function y(r) = g{r)e^^ r ^ at contact r = 1 
is found to be discontinuous as 

«/ SPS (i-) = ~— l -rrr (Bi2) 

lr a + 1/C 

p y SPS (i + ) = P y SPS (^) + -^j^. (B13) 

(" + 1/0 

Note that Eq. (B12) is the sticky limit of the PSW value py(l), Eq. (A4), [recall that y(r) is continuous at r = 1 
within the PSW] and is also the sticky limit of the PSW value py(X~), Eq. (A5). On the other hand. Eq. (B13) 
is the sticky limit of the PSW value y(X + ), Eq. (A5). Therefore, the discontinuity of y SFS (r) at r = 1 is a direct 
consequence of the discontinuity of the PSW cavity function at r = A. Both discontinuities are artifacts of the LPA. 
Again, this can be amended by an improved mLPA approach which is discussed in Sec. VII. 

APPENDIX C: LOW-DENSITY EXPANSION OF THE LPA 

Let us compare the LPA to order p with the exact results. From Eqs. (10)— (12) we easily get 

C = (oP + CiP 2 + 0(p 3 ) (CI) 

with 

Co = 7r, Ci = 7r(l-7A). (C2) 
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Upon inserting the result into Eqs. (14) and (15), and after some algebra, we find the correct zeroth order term 
9o( r ) = ffo xact ( r ) as gi ven m Eq. (17), and 

(l- 7r ) [l- 7 ^A-(r-l)l 



(1 + 77,-) l- 7 l^A + (r-l) y r - 7 2 7> 



5i (r) = 7r < 



(C3) 



2 - 27A - r, 
7(2 + 7) (r-2)-2 7 A, 
(2 + 2A-r) 7 2 , 
0. 

Comparison between Eqs. (C3) and (18) shows that the LPA reproduces the exact result for r > 1 + A. On the 
other hand, it fails to do so within the potential range. The differences between the first-order exact and LPA cavity 
functions are given by Eq. (19). 



< r < 1, 

1 < r < 1 + A, 

1 + A < r < 2, 

2 < r < 2 + A, 
2 + A<r<2 + 2A, 
2 + 2A < r. 
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Figure captions 

FIG. 1: (Color online) The PSW potential (top panel). The middle and bottom panels sketch the different behavior of the SW 
and PSW models, respectively. In the SW case there exists a hard core (black inner sphere) and an interaction range (light blue 
outer sphere) so two spheres on a line can either non interact (A) or attract each other as the corresponding interaction spheres 
overlap (B). As a consequence, different spheres cannot interchange positions on a one-dimensional line and the problem is 
analytically solvable. In the PSW the core is soft (red inner sphere) and hence we can have in addition to configurations (A) 
and (B) identical to the SW case, also the case where the internal cores overlap such as (C) and (D). Different spheres can then 
interchange position and the problem is a many-body one. 

FIG. 2: (Color online) Schematic phase diagram in the pa-ksT / e a space showing the region where LPA can be considered as 
reliable. The curves correspond, from top to bottom, to the cases (e r /e a , A) = (5,0.5), (5, 1), and (2,0.5). The points below 
each curve represent states where the relative deviation between the three routes to the pressure is smaller than 5%. 

FIG. 3: (Color online) Results for the radial distribution function g(r) versus r/a with A/a — 0.5, fcsT/e a = 1, e r /e a = 5, 
and pa = 0.2. Predictions from the modified LPA given by Eq. (54) (long dashed line) are compared with both MC results 
(solid line) and PY and HNC integral equations (short dashed and dotted lines, respectively). In the inset a magnification of 
the r < a region is shown. 

FIG. 4: (Color online) Same as in Fig. 3 at the instability threshold e r /e a — 2. All other parameters are as in Fig. 3. 

FIG. 5: (Color online) Same as in Fig. 3 at a higher density pa — 0.8. All other parameters are as in Fig. 3. 

FIG. 6: (Color online) Plot of (3pa vs pa for A/ a — 0.5, e r /e a = 5, and ksT/ta = 1 (top panel) and kBT/e a = 5 (bottom 
panel). Different curves refer to different routes. The symbols denote MC simulation results. 



FIG. 7: (Color online) Same as in Fig. 6, except that e r /e a = 2. 
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FIG. 8: (Color online) Plot of the Fisher-Widom transition line in the pcr/e a vs fc.gT/e a plane (top panel) and in the pa vs 
ksT/ta plane (bottom panel). Here A/cr = 1 and e r /e a = (oo, 10, 5). Note that, except in the SW case (e r /e a = oo), £/<r ^ (5p. 
Note also that in these cases the three routes to the pressure are not distinguishable one from the other on the graph scale. 



FIG. 9: (Color online) Same as in Fig. 8, except that A/cr = 0.5. 
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